Accepted for publication in The Astrophysical Journal 

Preprint typeset using 1^1^^ style cmulatcapj v. ^jllj^^ 



OBSERVATIONS OF RAPID DISK-JET INTERACTION IN THE MICROQUASAR GRS 1915+105 
David M. Rothstein,^ Stephen S. Eikenberry,^ and Keith Matthews^ 

Accented for publication in The Astrophysical Journal 

ABSTRACT 

We present evidence that ~ 30 minute episodes of jet formation in the Galactic microquasar 
GRS 1915+105 may sometimes entirely be a superposition of smaller, faster phenomena. We base this 
conclusion on simultaneous X-ray and infrared observations in July 2002, using the Rossi X-ray Tim- 
ing Explorer and the Palomar 5 meter telescope. On two nights, we observed quasi-periodic infrared 
flares from GRS 1915+105, each accompanied by a set of fast oscillations in the X-ray light curve 
(indicating an interaction between the jet and accretion disk). In contrast to similar observations in 
1997, we find that the duration of each X-ray cycle matches the duration of its accompanying infrared 
flare, and we observed one instance in which an isolated X-ray oscillation occurred at the same time as 
a faint infrared "subflare" (of duration 150 seconds) superimposed on one of the main flares. From 
these data, we are able to conclude that each X-ray oscillation had an associated faint infrared flare and 
that these flares blend together to form, and entirely comprise, the ~ 30 minute events we observed. 
Part of the infrared emission in 1997 also appears to be due to superimposed small flares, but it was 
overshadowed by infrared-bright ejections associated with the appearance of a sharp "trigger" spike 
in each X-ray cycle that were not present in 2002. We also study the evolution of the X-ray spectrum 
and find significant differences in the high energy power law component, which was strongly variable 
in 1997 but not in 2002. Taken together, these observations reveal the diversity of ways in which the 
accretion disk and jet in black hole systems are capable of interacting and solidify the importance of 
the trigger spike for large ejections to occur on ^ 30 minute timescales in GRS 1915+105. 
Subject headings: accretion, accretion disks — black hole physics — infrared: stars — stars: individual 
(GRS 1915+105) — X-rays: binaries 



1. INTRODUCTION 

GRS 1915+105 is one of the most fascinating objects in 
astrophysics today. It was one of the first systems in the 
Galaxy to be identified as a "microquasar" — an X-ray 
binary with relativistic jets which mimics some of the be- 
havior of quasars on a s maller and closer scale. (For a re- 
view of these objects, see lMirabel fc R.odr fguez"1999'and 
iFender 2004: f or a re view focusing on GRS 1915+105, 
see iFender fc B elloni 2 0041 ') Because timescales in mi- 
croquasars are expected to be a factor of ~ 10® shorter 
than in quasars (owing to the much smaller mass of the 
compact object which powers the system), microquasars 
are excellent laboratories for investigating accretion disk 
evolution and jet formation in black hole systems. 

GRS 1915+ 105 was first discovered a s a transient 
X-ray source (iCastro-Tirado et"al] Il994j) . and its ra- 
dio and infrared counterparts have since been identified 
((Mirabel et al. 1994). The system contains a ~ 14 M0 
black hole fed b y Roche lobe overflow from a K-M 
giant companion llGreiner. Cubv fc McCaughreaiJl2001t 
IGreiner et all I2001D . GRS 1915+105 displays extreme 
variability in many different wavebands and on many 
different timescales, making it unique among the X- 
ray binaries. Its uniqueness is likely due to its ex- 
tremely high accretion rate, which allows it to regularly 
reach luminosities unseen in other Galactic X-ray sources 
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Perhaps the most spectacular variability observed from 
GRS 1915+105 is seen in high resolution radio maps, 
which occasionally reveal the presence of resolved, bipo- 
lar ejections with flux densities up to ~ 600 mJy that 
form on timescales of weeks and m ove away from the sys- 
tem a t relativistic speeds (> 0.9c; Mirabel fc Rpdrigued 
IT99I lR,odrfguez fc Mirabell IT999t IFender et all Il999ar 
These events are referr ed to as "class A" ejections by 
lEikenberrv et al.1 1)2000|) to distinguish them from smaller 
"class B" radio and infrared flares on ~ 30 minute 
timescales which are also thought to correspond to jet 
ejection. In addition to the radio and infrared flares, 
GRS 1915+105 displays a broad range of variability in 
the X-rays, where it switches between maiiy diff erent 
states (reviewed extensively bv iBelloni et al.ll2000D and 
where quasi-periodic oscillations (QPO s) on timescales 
as fast as ~ 168 Hz have been observed (|R.emillard et alJ 
12003; McClintock & Remillard 2004). 

Multiwavelength observat ions of "class B " event s 
in GRS 1915+105 by iPoolev fc Fellde^ lll99l 
lEikenberrv et all l|1998albD . and iMirabel et alJ l|1998fl 
were the first to reveal the intimate link between accre- 
tion disk evolution and relativistic jet f ormation on short 
timescales in any black hole system. lEikenberrv et all 
|r998a|) observed a one-to-one correspondence between 
repeating X-ray variability cycles and infrared flares 
100 mJy dereddened) on timescales of ~ 30 minutes. 
The observations are consistent with a picture in which 
emptying and refilling of the X-ray emitting inner 
disk (e.g. Belloni et al., 1997a.h' ) coincides with the 
ejection of material into a jet, which radiates through 
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synchrotron emission to pro duce a flare ijFender et al.l 
119971 iPoolev fc Fended ITqQTI) . The infrared and X-ray 
hght curves decouple as the ejecta becomes causall y 
separated from the inner disk ijEikenberrv et alJll998a|) . 
and as the plasma moves out of the accretion disk 
plane it can rad iatively pump ernission l ines originating 
within the disk ijEikenberrv et al1ll998b(l . Each infrared 
flare appears to be accompanied by a delayed radio 
flare, perhaps indicating adiabati c expansion of the 
ejected cloud fMirabel ct al. 1998) or motion along a 
conical jet ifFtender fc Po oley 1998), while the jet itself is 
resolved by the Very Long Baseline Array as an optically 
thick synchrotron source of length ~ 20 AU that is 
variable on similar timescales as the X -rays and infrared 
ijDhawan. Mirabel fc Rodrfg uez"2000^. 

Though subsequent observations have revealed smaller, 
more complex infrared flari ng behavior in this source 
(e.g. lEikenberrv et al.1 12000(1 . the above represents the 
basic picture for the ejection of "class B" jets in 
GRS 1915-1-105. In this paper, we present multiwave- 
length observations of GRS 1915-1-105 on four nights dur- 
ing July 2002, when the source was undergoing another 
period of "class B" ejection. The similarities and dif- 
ferences between these observations and those obtained 
previously allow us to begin to map out how the many 
different X-ray states of GRS 1915-f 105 affect jet pro- 
duction in this unusual source. We conclude that the 
observed ejections may sometimes entirely be a superpo- 
sition of smaller, more complex phenomena, and that the 
appearance of a "trigger" spike in the X-ray light curve 
(seen in the observations of Eikcnbcrrv ct al. 1998a| but 
not in our July 2002 observations) is a key ingredient for 
large, infrared-bright ejections to occur on ~ 30 minute 
timescales. 

2. OBSERVATIONS 

We obtained infrared observations of GRS 1915-1-105 
on the nights of 2002 July 27-29 UT in the K (2.2/Ltm) 
band, using the D-78 camera at the Cassegrain focus of 
the Palomar Observatory 5 meter Hale telescope. X-ray 
observations (with coverage between ~ 2-100 keV) were 
obtained on July 27-28 and July 30 using the Propor- 
tional Counter Array (PGA) on the Rossi X-Ray Tim- 
ing Explorer (RXTE); fu rther details regarding this in- 
strument can be found in lGreiner. Morgan fc RemillardI 
(^96) and references therein. 

We configured the infrared camera to take 128 x 128 
pixel (16" X 16") images at a rate of one frame per sec- 
ond, with absolute timing provided by a WWVB radio 
signal from the National Institute of Standards and Tech- 
nology (~ 1 ms accuracy). We observed GRS 1915-1-105 
in this mode for ^ 8 hours on each of the first two nights 
and for ~ 1 hour on July 29. We processed each image 
by subtracting an averaged sky frame, dividing by a flat 
field, interpolating over bad pixels and filtering it in the 
Fourier domain to remove electronic pattern noise that 
corrupted many of the images. On July 27, we also ob- 
served the faint HS T/NIGMOS standa rd GSPG P182-E 
(star no. 9177 from lPersson et anil998D for absolute fiux 
calibration. 

A typical processed image is shown in Figure Q The 
field of view is large enough to capture GRS 19154-105 
and several nearby stars. We performed differential pho- 
tometry on GRS 19154-105 and the surrounding field 



Star A 



Reference k| 
stars ~ 



GRS 1915+105 



Fig. 1. — A typical 2.2fj,m infrared image of the field of 
GRS 19154-105 at 1-second time resolution, from our observations 
on 2002 July 27. North is up, east is to the left, and the field of 
view is 16" X 16". The two reference stars were used to correct the 
GRS 1915+105 and Star A light curves for atmospheric variability 
and changes in the point spread function. 



stars, measuring their fiuxes within a 1" radius soft- 
ware aperture. We used the two brightest field stars 
as a reference, averaging their light curves together and 
dividing the normalized result into the GRS 19154-105 
light curve to correct for atmospheric variability and 
changes in the point spread function. This procedure 
was also applied to "Star A" in Figure its flux den- 
sity remained steady at a value of 3.1 i nJy (consistent 
with the reported value in lFender et alj [r997'l and was 
used to determine the flux density of GRS 19154-105. 
The GRS 19154-105 light curve was then dereddened by 
Ak = 3.3 mag nitudes to correct for absorption in the 
Galactic plane l|Fender et alJll997|) . We chose Ak = 3.3 
to be consistent with the previous literature, although 
more recent work suggests that Ak = 2.2 is a better 
estimate for this source ( Fu chs. Mirabel fc Glared 1200^ 
IChapuis fc Cor bel 2004). If Ak = 2.2 is adopted, then 
all the infrared flux densities in this paper should be re- 
duced by a factor of ~ 3. 

For the X-ray analysis, we extracted light curves (at 
1-second time resolution) from PGA Standard-1 data 
using FTOOLS v5.2. We also extracted X-ray spec- 
tra in the ^ 3-25 keV range, using data from PGA 
Binned mode B_8ms_16A_0_35_H_4P and Event mode 
E_16us_16B_36_ls. We used standard procedures for re- 
sponse matrix generation, background estimation and 
subtraction, and correction for PGA deadtime. We 
then used XSPEG vll.2 to fit each spectrum with a 
standard model for black hole candidates consisting 
of a "soft" component (which peaks in the low en- 
ergy X-ray s) modeled as a multitemperature disk black- 
body (e.g. 'Mitsuda et al. 1984) and a "hard" compo- 
nent (which extends to the higher energy X-rays) mod- 
eled as a power law, both modified by hydrogen absorp- 
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Fig. 2. — Simultaneous infrared and X-ray light curves of GRS 1915-1-105 on (a) 2002 July 27 and (b) 2002 July 28. Both light curves are 
at 1-second resolution, but the infrared data have been smoothed to 5-second resolution in this plot. Typical 1-second ± la uncertainties 
are shown for the infrared; the X-ray uncertainties are assumed to be Poissonian and are too small to be seen on this scale. The infrared 
data were taken in the K (2.2/jm) band and have been dereddened by 3.3 magnitudes, while the X-ray data are in the ~ 2-100 keV energy 
range and are plotted in units of counts per second per Proportional Counter Unit (PCU) of the RXTE PCA instrument. Gaps in the 
X-ray light curves represent times when RXTE was not observing the source. The "subflare" and isolated "oscillation" in (a) are discussed 
in El 



tion fixed to a column density of 6 x 10^^ atoms cm ^ 
uno. Morgan fc Rem illard* 1999*). A systematic error 
of 1% was added to each spectrum before the fit was 
performed. 

3. LIGHT CURVES 

The resulting light curves from July 27-28 are shown 
in Figure 121 On each night, we detected quasi-periodic 
infrared flares with dereddened amplitudes of 30 mJy. 
The flares have durations of ~ 15-30 minutes and repeat 
on timescales of ~ 30-60 minutes. The accompanying 
X-ray light curves show GRS 1915-1-105 undergoing a se- 
ries of long "dips" and fast oscillations. The particular 
pattern of X-ray variability seen here has been observed 
many times in GRS 1915-1-105; it is the "class a" state 
defined by BcUoni ct al. (2000), characterized by oscilla- 
tions which grow longer and fainter with time before the 
source enters an X-ray dip of duration ~ 1000 seconds. 

Our observations on July 29 were limited to ~ 1 hour 



but showed two infrared flares similar to those seen on 
July 27-28. This suggests that the state we observed 
lasted for more than two days. The light curves from 
July 29 are not shown here because there were no ac- 
companying X-ray observations, but they are available 
upon request. X-ray observations on July 30 indicated 
that GRS 1915+105 was no longer undergoing oscilla- 
tions and had entered a perio d of hard, steady X-ray 
emission (the "class x" state of lBelloni erani200nn . 

3.1. Multiwavelength Features of the Light Curves 

As can be seen from Figure [2 each time the X-ray 
light curve of GRS 1915+105 switches from a dip into 
a period of oscillation, it is accompanied by an infrared 
flare. Furthermore, the flare appears to be triggered by 
a sharp X-ray "spike" of duration several seconds which 
is present during each transition. 

These features are broadly consistent with previous 
multiwavelength observations of GRS 1915+105 dur- 
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Fig. 3. — Simultaneous infrared and X-ray light curves of 
GRS 1915+105 at 1-secon d resolution during Augu st 1997, from 
data originally presented in lEikenberrv et al.l I1998ai) . The infrared 
data have been dereddened by 3.3 magnitudes before plotting. In 
(a), the infrared flare appears to end and the flux density decays 
back to its quiescent level while the X-rays are still undergoing os- 
cillations. In (b), there is an infrared flare beginning at ^ 8,000 
seconds which appears to decouple from the X-rays (though there 
is no X-ray coverage for the initial flare rise), but there is also 
a ~ 20 mjy "infrared excess" that continues until ~ 11,000 sec- 
onds and that is associated with an unusually long series of X-ray 
oscillations. 



ing "class B" ejections (e.g. lEikenberrv et al.1 Il998at 
^irabcl et al. 199 Si); in fa c t, man y of the X-ray states 
classified by .Belloni et al. ( 2000) in which the source 
transitions between a long dip and a period of oscilla- 
tions are now known to t rigger radio (and theref ore pre- 
sumably infrared) flares ijKlein- Wolt et alJl200l . 

What makes the current observations interesting is 
their high time resolution and broad coverage of many 
individual episodes of jet productio n, matc hed only by 
the similar observations of Eikcnberr v et alJ (ji998a). By 
looking at how the subtle differences in X-ray evolution 
between these two sets of observations lead to differences 
in the infrared behavior, we can begin to unravel the de- 
tails of the complicated evolution of GRS 1915-1-105. 

3.2. Origin of the Infrared Flares 

Infrared flares similar to those presented here have 
been observed many times from GRS 1915+105 and have 
been consis tently interpreted as synchrotron emission 
from a jet iFender et all Il997t iFender fc Poolevi Il998t 
Eikeriberrv et 998ab IMirahel et aLIII 998HTTeda et a,].l 
20021) . We propose that the flares in Figure 12 have the 
same origin. 

This interpretation is supported by 15 GHz observa- 
tions at the Ryle Telescope which were obtained sev- 
eral hours before our July 28 observations (G. Pooley 
2002, private communication). Radio flares with simi- 
lar timescales as the infrared flares and amplitudes of 
~ 10-20 mJy (above a baseline flux density of ~ 10 mJy) 
were observed. If the infrared flares we observed had sim- 



ilar radio counterparts, then the relatively flat infrared- 
to-radio spectrum (as well as the ~ 10^ to 10^° K bright- 
ness temperatures of the radio flares) indicates a nonther- 
mal origin. 

4. INFRARED /X-RAY CORRELATION 

The most interesting result from these observations 
may be the strong correlation between the duration of 
each infrared flare and the duration of its accompany- 
ing X-ray oscillation period. All the flares have similar 
rise times (~ 200-300 seconds), but each flare appears 
to "wait" to return to its quiescent level until the X-rays 
stop oscillating. This gives rise to asymmetric profiles 
for some of the longer flares, particularly those on July 
27. 

How does this behavior compare to that seen in previ- 
ous observations of GR S 1915-1-105? The fla res observed 
in August 1997 bv .Eikenberrv et all l)1998aj) were a few 
times stronger than those presented here, with peak am- 
plitudes as high as ~ 200 mJy dereddened. A correlation 
between the infrared flare duration and the X-ray oscil- 
lation period was not directly evident in 1997, and a few 
flares appear to have decayed almost completely while 
the X-rays were still undergoing strong oscillations. An 
example of this behavior is shown in Figure 

However, there is some evidence for a smaller contribu- 
tion to the infrared light curve (on the ~ 20 mJy level) 
associated with the 1997 X-ray oscillations. In particu- 
lar, a long series of oscillations occurred on August 15 at 
the same time as a period of "infrared excess" following a 
flare, in which the GRS 1915-f 105 flux density remained 
steady at a level above it s quiescent value (see Figure|3|D). 
lEikenberrv et all 1)20001) showed that if each X-ray oscil- 
lation in this series had a '--^ 5-10 mJy infrared flare asso- 
ciated with it, then the superposition of these flares could 
reproduce the observed infrared excess. Their conclusion 
was based on the spectral similarity of the 1997 oscilla- 
tions to isolated X-ray oscillations in July 1998, each of 
which had an accompanying, time-resolved ~ 5-10 mJy 
infrared flare (the "class C" flares). 

The observations presented in o ur current work 
strengthen the argument made by lEikenberrv et aP 
(2000). The key piece of evidence comes from the third 
infrared flare on July 27 (Figure |2ti) . This flare contains 
clear evidence for a "subflare" — a secondary peak that 
occurs after the flare has begun to decay. Furthermore, 
the subflare is associated with an isolated oscillation in 
the X-ray light curve that occurs after an unusually long 
delay between it and the previous set of oscillations. It is 
therefore plausible that each X-ray oscillation has a faint 
infrared subflare associated with it, but in general the os- 
cillations are so closely spaced in time that the subflares 
blend together to create the illusion of one continuous 
flare.''' 

Unlike in the 1997 data, where these faint flares sim- 
ply contribute an "infrared excess" on top of a larger 
episode of flaring behavior, we propose that the ^ 30 
minute flares in 2002 are composed entirely of these su- 
perimposed events. This would immediately explain the 
correlation between the infrared flare duration and the 

* Note that there are a few other infrared flares in Figure |5| 
that appear to have very faint subflares, but none of these events 
have good enough X-ray coverage to determine whether they are 
associated with isolated X-ray oscillations. 
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Fig. 4. — The number of X-ray oscillations in each rsj 30 minute 
cycle plotted against the integrated flux density of the infrared 
flare accompanying that cycle, from our observations on 2002 July 
27—28. Data from all eight infrared flares that had uninterrupted 
X-ray coverage are shown in this flgure, plotted as points. (Note 
that two of the points with 12 X-ray oscillations are difficult to 
discern because they almost completely overlap.) The dotted line 
shows the best linear fit to the data, and the dashed line shows the 
best linear fit constrained to pass through the origin. 



X-ray oscillation period seen in 2002; the flare is only 
sustained by the continual production of subflares asso- 
ciated with the X-ray oscillations. 

To test whether this idea is feasible, we show in Figure 
01 the number of X-ray oscillations in each ~ 30 minute 
cycle plotted against the integrated flux density of the 
infrared flare accompanying that cycle. (The integrated 
flux densities were calculated by smoothing the light 
curve to 30 second resolution, determining the start and 
end times of each flare by eye, subtracting out the base- 
line and integrating the result numerically.) It is clear 
from Figure 01 that X-ray cycles with more oscillations 
produce more infrared emission, consistent with our sug- 
gestion. Figure 01 also shows linear fits to the data; the 
best fit line has a slope of ~ 1.9 Jy sec oscillation"^ and 
a vertical intercept of ~ 3.6 Jy sec, while the best fit 
line constrained to pass through the origin has a slope of 
^ 2.2 Jy sec oscillation"^. Both lines appear to be rea- 
sonable fits, and if we assume a plausible value of ~ 1.5 
Jy sec for our measurement errors in determining the 
integrated flux densities, then both fits are found to be 
statistically acceptable, with reduced chi-squared values 
of xt 0.4 and 1.1, respectively. 

We also tried fitting a power law to the data in Fig- 
ure 01 and found a slightly nonlinear relationship (power 
law index 0.8), but the quality of the fit was not sig- 
nificantly better than that of the linear fit {xt ~ 0.34). 
Therefore, the simplest conclusion we can draw from our 
observations is that there is a roughly linear relationship 
between the number of X-ray oscillations in each ^ 30 
minute cycle and the strength of the infrared flare ac- 
companying that cycle. In addition, we tested whether 
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Fig. 5. — Simultaneous infrared and X-ray light curves of 
GRS 1915-1-105 at the time of the sixth infrared flare observed 
on 2002 July 27 (Figurel^i). The infrared data (dereddened by 3.3 
magnitudes) are plotted as points at 1-second resolution; the curve 
in the top panel shows the results of our simulation (discussed in 
CT . in which the infrared flare is modeled as a superposition of 
small flares, one for each X-ray oscillation. 



the infrared flare strength might actually be related to 
the total amount of X-ray emission during the oscilla- 
tions (rather than the number of oscillations) and found 
no evidence to support this; if the infrared flare strength 
is plotted against the integrated X-ray flux in each cy- 
cle (above a baseline level of 1800 counts/sec/PCU), the 
correlation has a larger scatter than that shown in Fig- 
ure 01 and both linear and power law fits to the data are 
poor ixl>2). ^ 

The simplest interpretation of the linear relationship 
in Figure 01 is that each X-ray oscillation contributes a 
constant amount of infrared emission to the overall light 
curve, in the form of a faint flare. In addition, the fact 
that the vertical intercept of our fit is small (equivalent 
to an infrared excess of only a few mJy over the typi- 
cal duration of a flare) and possibly consistent with zero 
supports our suggestion that the infrared events are com- 
posed entirely of these superimposed faint flares. Fur- 
thermore, the slope that we measure from our fit indi- 
cates that each X-ray oscillation should be responsible for 
~ 2 Jy sec of infrared emission. This number is slightly 
larger than, but still consistent with, the integrated flux 
density in the sub flare we observed (between ~ 0.5 and 
2 Jy sec, depending on how we choose to distinguish the 
"subflare" from the underlying main flare). It may also 
be reasonable to assume that an average X-ray oscilla- 
tion triggers a slightly larger subflare than the one we 
observed, since this subflare was associated with one of 
the weakest oscillations in its series. 

To further test our idea, we attempted to simulate the 
infrared light curve as a superposition of faint subflares. 
We began our simulation with a constant infrared back- 
ground level of ~ 70 mJy and added one faint flare to 
this background for each oscillation in the X-ray data. 
We modeled each flare as a Gaussian with fixed ampli- 
tude and full width at half maximum (FWHM). Based on 
the observed properties of the July 27 subflare, we chose 
the position of each Gaussian so that its peak occurred 60 
seconds after the peak of its accompanying X-ray oscil- 
lation. We performed several simulations with different, 
fixed values of the amplitude and FWHM, since these 
quantities were hard to determine from the observed sub- 
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Fig. 6. — Comparison of one cycle of jet formation in GRS 1915+105 from observations in (a) 1997 and (b) 2002. The top two panels in 
each column show the infrared (dereddened by 3.3 magnitudes) and X-ray light curves, respectively (at 1-second resolution, but the infrared 
data have been smoothed to 5-second resolution in this plot); the dotted horizontal line in (b) shows the baseline infrared level. The bottom 
four panels show the results of our X-ray spectral fitting at 4-second resolution, plotted along with ± la uncertainties; these include the 
unabsorbed 2—25 keV flux from the multitemperature disk blackbody and power law components of the spectrum, the temperature at the 
inner edge of the accretion disk, the radius of the inner edge of the accretion disk, and the power law index. The vertical shaded regions 
numbered 1 through 7 show the parts of the light curve from which the spectra plotted in Figure |S] were extracted. 



flare. 

In general, these simulations did not fit the data 
well; the superimposed flares we produced reached peaks 
which were initially too high or dropped to values which 
were too low by the end of the ^ 30 minute cycle. We 



therefore decided to allow the FWHM of the flares to 
vary based on the duration of each X-ray oscillation. We 
chose a simple linear scaling, in which the FWHM of each 
flare was taken to be proportional to the rise time of its 
accompanying X-ray oscillation. Results of this Simula- 
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tion for the sixth infrared flare on July 27 are shown in 
Figure El in this case, we fixed the amphtude of each 
Gaussian to be 8 mJy and set each Gaussian's FWHM 
to be a factor of 5 greater than the rise time of its accom- 
panying X-ray oscillation. This choice guaranteed that a 
typical oscillation in each X-ray series, with a rise time 
of ~ 50 seconds, had an associated infrared flare with a 
FWHM of 250 seconds and an integrated flux density of 
2 Jy sec, consistent with the slope of the fit in Figure 

Simulations such as the one in Figure[Sldo a good job of 
reproducing the observations, especially given the gross 
approximations we are making. Even so, some difficulties 
are encountered. The July 27 subfiare is not reproduced 
using the above procedure; its accompanying X-ray os- 
cillation has an extremely long rise time, and therefore 
the simulation tries to model it with an unrealistically 
large FWHM. In general, the simulations work better 
at the beginning of each ^ 30 minute cycle than at the 
end, suggesting that there is residual infrared emission 
associated with the end of each X-ray cycle that is not 
captured in our model. Nonetheless, the results shown in 
Figure are fairly typical, and we can be confident that 
the 2002 fiares consist primarily of superimposed faint 
flares associated with each X-ray oscillation. 



5. X-RAY SPECTRAL EVOLUTION 

We can compare our X-ray light c urves and spectra 
to thos e obtained in August 1997 by lEikenberrv et al.l 
()1998aj) . Though both sets of observations show tran- 
sitions between X-ray dips and oscillations, the details 
of the beh avior are different; in the classiflcation of 
iBelloni et al. (2000), GRS 1915+105 was in the "class 
/3" state in 1997 and the "class a" state in 2002. 

In Figure ini we present one cycle of jet formation from 
each of the 1997 and 2002 observations. Several pa- 
rameters obtained from our spectral fltting are plotted 
at 4-second resolution: the unabsorbed 2-25 keV flux 
from the multitemperature disk blackbody and power 
law components of the spectrum, the temperature at the 
inner edge of the accretion disk Tin, the radius of the 
inner edge Rm, and the power law index 7 (where the 
flux density due to the power law, at photon energy E, 
is oc E~^). The fits shown in Figure have typical re- 
duced chi-squared values of ^ 1-2, and all fits with 
> 2, as well as those for which XSPEC was unable to 
return uncertainties on the fit parameters, were excluded 
from the plot. Note that Rm is calculated from the disk 
blackbody normalization {Rm (x \fN , where N is the 
normalization) assuming a distance to GRS 1915+105 of 
11 kp c and an accretion d isk inclined 66° to the line of 
sight l|Fender et alJll999aj) . but the values derived from 
this procedure should not be viewed as accurate esti- 
mates of the physical size o f the disk (see ti5.2.2l and 
iMerloni. Fabian fc Rossll2000|) . 

^ Note that by allowing the FWHM to vary, we are also varying 
the amount of infrared emission associated with each X-ray oscil- 
lation. However, this does not contradict the linear relationship 
found in Figure |1] it turns out that each X-ray cycle has about 
an equal number of oscillations with rise times above and below 
the typical value of ~ 50 seconds, so the stronger infrared flares 
may "cancel out" the weaker ones and still lead to a roughly linear 
relationship. 



5.1. 1997 and 2002 State Changes 

Figure ini shows that in both 1997 and 2002, the X-ray 
spike which appears near the beginning of the infrared 
fiare (at ~ 2, 200 seconds in 1997 and ~ 24, 000 seconds 
in 2002) represents a "state change" within the accretion 
disk, in the sense that it initiates changes that persist af- 
ter the spike is finished. During the X-ray dip, the source 
is in a spectrally hard state in which the disk blackbody 
component is nearly absent, but the increase in Tin and 
in the disk emission which occur at the spike indicate 
that the disk becomes hot and visible at this point and 
remains that way throughout the X-ray oscillations. 

In 1997, there is a prolonged period of a few hundred 
seconds after the X-ray spike when the inner accretion 
disk has "turned on" but before the period of fast oscilla- 
tions. It is clear that the infrared flare begins well before 
the oscillations, so although we may be able to ascribe 
some of the infrared emission to blended together "sub- 
flares" associated with each oscillation (as discussed in 
21 , there must be some other mechanism at work which 
leads to the initial rise of the infrared light curve right 
around the moment of the spike and allows it to reach a 
peak amplitude of ~ 100 mJy. 

In 2002, however, the spike seems to represent the first 
oscillation in a series; it is spectrally similar to the os- 
cillations which immediately follow it, and differs mainly 
in that it is the strongest and fastest oscillation in the 
group. This is consistent with the idea that the flares ob- 
served in 2002 consist entirely of superimposed ^ 5-10 
mJy subflares associated with each oscillation, and that 
no abnormally large ejection occurs at the moment of the 
flrst spike. 

The differences between the 1997 and 2002 behavior 
can be more clearly seen in Figure where we have 
plotted the results of X-ray spectral fitting performed 
at high time resolution (1 second). Again, all fits with 
> 2, as well as those for which XSPEC was unable to 
return uncertainties on the fit parameters, were excluded 
from the plot. Figure da, shows the initial X-ray spike in 
1997, Figure[7lD shows the initial X-ray spike in 2002 and 
the oscillations which immediately follow it, and Figure 
Et shows a typical set of oscillations in 1997. 

It is clear from Figure [3 that although the details of 
the 1997 and 2002 oscillations are different, there are 
also many similarities: both consist primarily of peri- 
odic changes in the disk emission and Tin- The initial 
oscillation in 2002, which signals the beginning of the 
disk activity and coincides with the rise of the infrared 
flare, is no different in this regard. Figure however, 
shows that the initial spike in 1997 (which we refer to as 
the "trigger" spike) is an entirely different phenomenon. 
During this spike, Ti„ rises monotonically to ^ 1.8 keV, 
where it remains, while the power law component of the 
spectrum does not vary noticeably during the spike but 
changes sha rply as the spike ends. Sim ilar behavior was 
observed bv lMjidiari fc Bellonil l|200riil . but our spectral 
fits are performed at higher time resolution and therefore 
confirm that the power law does not change significantly 
even during the last few seconds of the spike's rise. 

5.2. 1997 Trigger Spike 

Given the importance of the 1997 trigger spike and its 
possible connection to large ejections in GRS 1915+105, 
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Fig. 7. — Comparison of the X-ray spectral evolution of GRS 1915+105 during (a) the 1997 trigger spike (from FigurelGkl. (b) the initial 
spike in 2002 and the oscillations which immediately follow it (Figure 03), and (c) a typical set of oscillations in 1997 (Figure |^). The 
top panel in each column shows the X-ray light curve (at 1-second resolution), and the bottom four panels show the same X-ray spectral 
fitting parameters plotted in Figure |H] but the fits shown here were performed at 1-second resolution. It is clear from this figure that the 
trigger spike in (a) is spectrally different from the other oscillations in 1997 and 2002. 



it is useful to examine its evolution in a model- 
independent way. This is shown in Figure IHt, where 
we plot three X-ray spectra associated with the trigger 
spike (the parts of the X-ray light curve from which these 
spectra have been extracted are indicated by the vertical 
shaded regions numbered 1 through 3 in Figure IHt^). We 
extracted the spectra at a time resolution of 16 seconds 
using PCA Staiidard-2 data, which have higher energy 
resolution than the Binned and Event mode data used 



for the fitting procedure. 

As can be seen from Figure between regions 1 and 
2 (as the spike begins) there is virtually no change in 
the hard (> 20 keV) X-ray emission, but the soft X-rays 
brighten dramatically. After the drop at the end of the 
spike, however (region 3), the spectrum steepens sharply, 
and the hard X-rays become nearly undetectable. The 
soft X-rays in region 3 retain a similar spectral shape as 
in region 2, but decrease in intensity. 
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In this model-independent description, the trigger 
spike consists of: 

1 . A factor of a few increase in the soft X-ray emission 
on a timescale of ~ 20 seconds or longer. 

2. A dramatic, factor of a few decrease in the hard X- 
ray emission and less dramatic decrease in the soft 
X-ray emission, on a timescale of a few seconds. 

How do these changes in the soft and hard X-ray emis- 
sion during the trigger spike relate to each other phys- 
ically? There are several possible ways to understand 
this, depending on how much we believe the accuracy 
of the multitemperature disk blackbody model that we 
have employed above. 

5.2.1. Nonthermal Models for the Spike 

Although our spectral fitting indicates that the soft X- 
ray excess during the trigger spike consists of thermal 
emission from the accretion disk, this is not an ironclad 
conclusion. The energy resolution of our data is quite 
poor, and it is possible that other models could repro- 
duce this brief burst of soft X-ray emission. Before the 
spike, the accretion disk component of the spectrum is 
extremely weak, and in fact, we are able to obtain accept- 
able fits at 1-second resolution {xt ~ 1-4) for a model 
that only includes a power law and no accretion disk. 
Thus, it is possible to conclude that the spike is sim- 
ply a "blip" on an otherwise smooth rise of the disk flux 
that corresponds (along with the smooth rise in Ti„) to 
the accretion disk "turning on." In this picture, when 
the 1997 X-ray light curve transitions from a dip into 
a period of oscillations, the thermal emission from the 
disk gradually becomes visible (as it does in 2002). At 
some point during this transition, however, a ~ 10 sec- 
ond nonthermal burst of soft X-rays is released, creating 
the spike. This interpretation certainly seems plausible 
in light of the second panel of FigureC^, which shows the 
disk flux rising smoothly before the spike and finishing its 
smooth rise afterwards; the spike can easily be imagined 
as a sharp feature superimposed on an otherwise smooth 
event. 

If this is the case, then when the nonthermal soft X-ray 
emission disappears at the end of the spike, the nonther- 
mal hard X-ray emission (from the power law compo- 
nent of the spectrum) is also observed to drop (Figure 
|S^). The origin of the power law in X-ray binaries is a 
subject of some debate; it is generally assumed to rep- 
resent either inverse Compton emissio n from a corona of 
relativistic electrons fPout anen'" 1 998*) which may consti- 
tute the base of a jet (Fcndei^et al,, 1999b), or the high 
frequency tail of syn chrotro n ernission from the iet i tself 
ijMarkoff. Falcke fc Fender. 200lllMa7koff et al.l200^ . In 
either case, it comes from a region close to the accre- 
tion disk, so the decrease in the power law emission at 
the end of the spike indicates that changes are occurring 
immediately outside of the accretion disk at this time. 
Meanwhile, the infrared flare also begins, as energy is 
dissipated in the jet further downstream from the disk. 

To test the nonthermal scenario, we tried fitting the 
soft X-ray emission during the spike as a power law (that 
is, we fit two power laws to the spectrum, the original 
one which extended to the hard X-rays, and a stronger, 
steeper one which covered the soft X-ray excess) but were 
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Fig. 8.— X-ray spectra of GRS 1915+105 during (a) the 1997 
trigger spike in Figure ISh. (b) the 1997 oscillations in Figure 15k. 
and (c) the 2002 oscillations in Figure 15b. The spectra are num- 
bered 1 through 7, corresponding to the regions of the X-ray light 
curve from which they were extracted (see Figure I^J. The spectra 
were obtained at 16-second resolution and are plotted as points 
with ± 1(T uncertainties; the curves in (b) and (c) show the power 
law component of the fit to each spectrum, modified by hydro- 
gen absorption (this component accounts for all of the hard X-ray 
emission and a fraction of the soft X-ray emission). 



not able to obtain acceptable fits at 1-second resolution. 
This implies the need for a "curved" spectral component 
to explain the emission during the spike; however, there 
are many possible sources for such a component besides 
thermal emission from the disk. 

The above suggestions are tentative, and much more 
data would be required to confirm them. Neverthe- 
less, a logical picture can be constructed in which the 
spike does not originate from thermal emission in the 
accretion disk. This picture would be consistent with 
several recent models, inc l uding the "magnetic bomb" 
fEikenbcrrv & van Puttcn 2003) and "magnetic flood" 
(,Tagger et a l. 20iI4) models for G RS 1915+105 and the 
more general model proposed by iLivio. Pringle fc Kind 
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1)2003(1 . These models all evoke the existence of large 
scale magnetic fields during the X-ray dip and associate 
the state transition at the end of the dip with the de- 
struction of these magnetic fields, perhaps through a re- 
connection event that produces a ~ 10 second burst of 
nonthermal emission along with an ejection of material. 

It is worth pointing out that our 2002 observations 
strengthen some of the above models, since they show 
that the large ejections 100 mJy infrared flares which 
decouple from the X-ray light curve) only seem to occur 
in the presence of the trigger spike. However, they also 
show that the ~ 30 minute state transitions can occur 
without an initial large infrared event. Why is this the 
case? If the spike and ejection represent the destruc- 
tion of the magnetic field which triggers the state tran- 
sition, then why do state transitions sometimes seem to 
require a trigger spike and ejection ~ 300 seconds before 
the fast oscillations, but other times occur without them 
and jump immediately into the oscillations? Ideally, a 
model for the GRS 1915+105 behavior should answer 
these questions. 

5.2.2. Thermal Accretion Disk Models for the Spike 

If we instead elect to believe that the soft X-ray excess 
during the trigger spike is due to thermal emission from 
the accretion disk, we need to explain how the disk in- 
teracts with the corona or jet (which is responsible for 
the power law component of the spectrum) in order to 
p roduce t he ob served emission. 

lYadavl 1(20011) suggested a general picture for the spike 
in which the soft X-ray luminosity increases until it 
reaches a critical value, thus triggering the ejection and 
subsequent changes in the light curve and spectrum. In 
this picture, all of the "action" takes place at the end 
of the spike; at this time, the X-ray power law and in- 
frared light curve change in a way that is consistent with 
jet ejection (just as they did in the nonthermal spike sce- 
nario of ^5.2.11) , but this change is now accompanied by a 
sharp drop in the accretion disk luminosity. This scenario 
has some advantages for observations in which the trigger 
spike does not have a well-defined beginning but rather 
seems to rise gradually to its maximum throughout the 
preceding X -ray dip (see, for e xample, the observation 
presented in lMidiari fc Bellonil l2003l . 

What could cause the accretion disk luminosity to de- 
crease at the end of the spike, at the same time as the 
power law component changes? One possibility is that 
the drop in disk luminosity indicates a sharp increase in 
the fraction of the accretion power which is being dis- 
sipated in the GRS 1915-1-105 jet. Although detailed 
modeling of the disk spectrum would be required to con- 
firm this suggestion, it is consistent with the changes ob- 
served in the disk parameters at the end of the spike (as 
discussed below), and it is also consistent with the ob- 
served changes in the X-ray power law and infrared light 
curve, which show that material may be ejected into the 
jet at the end of the spike. 

To study in detail the changes in the disk parameters at 
the end of the spike, it is useful to examine Figure[3i. As 
can be seen from this figure, it is only the fiux of radiation 
coming from the accretion disk which changes dramati- 
cally at the end of the spike, not the basic shape of the 
disk spectrum; the end of the spike occurs because of a 
sharp drop in i?i„ over a few seconds, while Tin increases 



over a slower timescale (~ 20 seconds) throughout the 
spike. 

This behavior is puzzling. In general, changes in the 
disk parameters occur on similar timescales, in which 
case an increase in Tin coupled with a decrease in Rin can 
be understood by a model in which the physical condi- 
tions in the outer part of the disk remain similar but the 
inner part of the disk has "turned on," either by refilling 
with matter o r by begi nning to emit detectable radiation 
((Belloni et a l."1997a bl. This scenario has been success- 
fully simulated by Watarai fc Mineshigg (2003); the de- 
crease in Rin is due to the fact tha t at high accretion 
rates (where the "slim disk" model of lAbramowicz et alJ 
119881 applies), significant thermal radiation can emerge 
from inside the last stable circular orbit around the black 
hole (Watarai et al. 2000). However, the fact that in our 
observations the drop in Rin occurs rapidly at the end 
of the spike, without any corresponding rapid changes 
in Tin, is more difficult to understand, and it is there- 
fore necessary to look for other interpretations for what 
is happening to the radiation from the accretion disk at 
t he end of the spik e. 

IRodriguez et alJ l(2002() point out that when the nor- 
malization of the multitemperature disk blackbody is ob- 
served to drop, the simplest interpretation is that the 
effective area of the emitting region of the disk has de- 
creased; the calculation of Rin from the normalization 
is model-dependent and not necessarily reliable. A drop 
in the fitted value of Rin can therefore be interpreted 
as the appearance of a hot spot or spiral shock pattern 
in the disk (Rodriguez et al. 2002), or more generally 
as a situation where some representative portion of the 
disk stops emitting radiation which we re ceive in the soft 
X-ray band (see als o Mun o et alJl200l1) . In this vein, 
Merloni. F abian fc Rossi (j200l ?) showed that if a more re- 
alistic model of the disk flow is used to simulate the X-ray 
spectrum, then changes in several parameters, including 
the fraction of the accretion power dissipated in regions 
outside of the disk (for example, in a corona or jet), could 
lead to apparent changes of Rin by a factor of ~ 2 when 
the multitemperature disk blackbody model is used to fit 
the spectrum. 

It is therefore plausible that the end of the spike co- 
incides with a transfer of accretion power from the disk 
(where it is dissipated thermally) to the jet (where it 
initiates changes in the X-ray power law and infrared 
flare) . During the ~ 300 second break between the spike 
and the subsequent oscillations, the disk then settles into 
the standard "slim disk" state, where the oscillations can 
occur. However, as in ti5.2.1l it is still not clear from this 
scenario what physical parameter triggers the sharp tran- 
sition and jet ejection in our 1997 observations but not 
in 2002. 

5.2.3. The Spike in the Context of Other X-ray Binaries 

iFender. Belloni fc Gallol l|2004D have proposed a uni- 
fled model for jet formation in X-ray binary systems. The 
model is "unified" in the sense that the trigger spike in 
GRS 1915-1-105 is shown to have counterparts in several 
other sources, where a radio flare is produced around 
the time of a peak in the X-ray light curve that oc- 
curs durin g the transition from a hard to soft spectral 
state. Fender. Belloni fc Gallol p004) proposed a phys- 
ical model in which a steady jet exists during the spec- 
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trally hard X-ray state preceding the spike and the speed 
of the jet increases sharply as the X-ray spectrum be- 
comes softer; the jet is then quenched completely once 
the softness reaches a critical value. Thus, synchrotron 
flares are produced via internal shocks in the preexist- 
ing jet around the time of the X-ray spike, when the jet 
velocity is increasing the fastest. 

Our observations are relevant to this model in several 
ways. First, we have shown that the sharp X-ray trigger 
spike observed in 1997 is associated with the appearance 
of large infrared flares that do not occur in 2002, when 
the spike is absent. This would seem to b e evidence in fa- 
vor of the lFender. Belloni fc Gallol lj2004j) model. On the 
other hand, a hard-to-soft state transition clearly occurs 
in 2002 as well, and the first oscillation in each 2002 cycle 
might therefore meet the definition of a trigger spike in 
this sense (see also ^7.1\\ . If so, the model would need to 
explain why the 1997 trigger spike is associated with a 
^ 100 mJy infrared flare, while the 2002 spike is only as- 
sociated with a ~ 5-10 mJy flare that is rivaled, at least 
in the infrared, by the flares that immediately follow it. 

There are several reasons why the strength of the flare 
might vary. As suggested bv lVadawale et all lj2003j) . the 
flare strength may be related to the total amount of ma- 
terial in the path of the faster jet — that is, if more 
material was ejected during the spectrally hard X-ray 
dip that preceded the spike, the resulting flare might be 
stronger. However, our observations do not show any 
particular relationship between the duration of each X- 
ray dip and the strength of the subsequent infrared flare. 
Therefore, there is no evidence to suggest that this effect 
is responsible for the stronger flares in 1997 as compared 
to 2002. 

Another mechanism for varying the str ength of the 
flare in the iFender. Belloni fc G& M ('2004') model is to 
change the speed of the jet. Render fc Belloni (200^ 
suggested that the proposed increase in jet speed as the 
X-ray spectrum softens may be due to the inner edge 
of the accretion disk moving closer to the black hole, 
where the escape velocity increases rapidly. If we take 
the results of our X-ray spectral fits in Figure [7| at face 
value, we can attempt to test this suggestion, by seeing 
whether there are significant differences in the evolution 
of Rin between 1997 and 2002. 

It can be seen from Figure [7\ that although i?i„ de- 
creases during the spike, it does so by a similar order 
of magnitude in 1997 as in 2002, suggesting that there 
may not be any large differences in the jet velocity. On 
the other hand, the value of Rm just before the spike 
appears to be slightly larger in 1997 than it is in 2002. 
Therefore, if the argument of lFender fc Be lloni (2004) is 
correct, the velocity of the steady jet is slower in 1997, 
leading to a greater velocity differential and stronger in- 
ternal shock when the highly relativistic jet is released 
during the spike. In addition, the decrease in i?i„ hap- 
pens on a faster timescale in 1997 than in 2002, which 
could mean that the jet velocity increases faster, leading 
to an internal shock in 1997 which occurs closer to the 
disk and therefore in the peak infrared-emitting region 
of the jet. 

However, this interpretation carries with it many 
caveats since, as we have pointed out in MS. 2. 21 the ob- 
served change in Rin may not correspond to an actual 
change in the inner radius of the disk, and even if it 



does, it only corresponds to the radius of the thermally 
emitting part of the disk (i.e. the region that is optically 
thick), which is not necessarily the same as the region 
from which the jet is ejected. Furthermore, there is a 
considerable range of X-ray behavior in both the 1997 
and 2002 observations, so any conclusions that apply to 
the data shown in Figure [7| do not necessarily apply as 
strongly to the "class a" and "class /9" states as a whole. 

An other way in which our observ ations are relevant to 
the Fe nder. Belloni fc Gallol (2004) model concerns the 
infrared subflare and our suggestion that each X-ray os- 
cillation in 2002 is associated with a similar amount of 
infrared emission (see 21 ■ If these faint infrared flares 
are interpreted as emission from the jet, then the jet in 
GRS 1915+105 may be active during the time of the 
soft X-ray oscillatio ns, in contrast to the suggestion of 
IFender. Belloni fc G allo ( 2001). On the other hand, we 
cannot conclusively prove that the subflares represent jet 
emission (see fj^Jl, and even if they do, it is not clear ex- 
actly where the "di viding line" between jet and non-jet 
states in the lFender. Belloni fc Gallo (2004) model actu- 
ally exists; it is possible that the low points of the 2002 
X-ray oscillations (where Tin ^1-5 keV) are spectrally 
hard enough to produce a steady jet. 

In general, our observations do not contradict the 
[Fender, Belloni fc Gallo ( 2004) model at this stage of its 
development; they will, however, likely be important in 
testi ng a more quantitative model based on the ideas 
that IFender. Belloni fc Gallcl (|2004|) proposed. Further- 
more, any successful model will need to explain what 
parameter triggers the sharp X-ray state transitions and 
infrared-bright ejections that are observed in 1997 but 
not in 2002. 

5.3. 1997 and 2002 Oscillations 

Perhaps the most interesting difference between the 
1997 and 2002 X-ray spectra is in the evolution of the 
power law component. In 1997, the power law index 
undergoes much more dramatic changes than in 2002, 
especially during the period of oscillations (see Figures IHl 
and[7|). If this difference is real, it would suggest that sig- 
nificant changes are occurring in the corona or jet during 
the 1997 oscillations that do not take place during 2002, 
which is somewhat puzzling since both sets of oscillations 
produce similar "class C" infrared flares (see 

We investigated the possibility that the observed varia- 
tions in the power law index were artificial. For example, 
if the multitemperature disk blackbody model provided 
a poor fit to the soft X-ray data, this could force the 
fitting routine to alter the power law component in or- 
der to fix the problems with the disk model on the soft 
X-ray end. Since the power law contributes significantly 
to both the soft and hard X-rays, the model spectrum 
might then fit the hard X-rays poorly, which would not 
necessarily be noticed in our fitting routine because the 
hard X-rays only provide a small percentage of the total 
fiux and contribute little to the chi-squared value of the 
fit. Therefore, it might be possible to obtain changes in 
the fitted power law index even if there were no changes 
in the actual flux of hard X-ray photons, which comprise 
the bulk of the power law shape. 

To get a handle on whether or not this is occurring, 
we show plots of X-ray spectra from the 1997 and 2002 
oscillations in Figure ISb and IS]:. The parts of the X-ray 
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Fig. 9. — Simultaneous infrared and X-ray light curves of GRS 1915+105 at the time of the infrared "subflare" observed on 2002 July 
27 (Figure|^). The X-ray data are at 1-second resolution, and the infrared data (dereddened by 3.3 magnitudes) have been smoothed to 
11-second resolution to increase the signal- to- noise ratio; they are plotted along with their it Icr uncertainties. The dotted vertical line 
marks the time of the X-ray spike, and the dashed horizontal line marks the median value of the infrared light curve calculated between 
6,030 and 6,110 seconds (the bold region of the dashed line). The subflare appears to rise above this value 30 seconds before the X-ray 
spike begins. 



light curve from which these spectra have been extracted 
are indicated by the vertical shaded regions numbered 
4 through 7 in Figure |H| As before, we extracted the 
spectra at a time resolution of 16 seconds using PCA 
Standard-2 data, in order to obtain the highest possible 
energy resolution. 

It is clear from Figure |H1 that the spectrum at energies 
> 20 keV does not change significantly during the 2002 
oscillations but does during 1997, where the hard X-ray 
flux appears to drop by a factor of a few as the light curve 
enters the low point of its oscillation (the low point of the 
oscillation is an ex ample of an X-ray "soft dip" discussed 
in e.g. Eikcnbcrrv et aljl20fl(lD . Furthermore, the solid 
lines in Figure |H1 show the power law fit to the data; at 
high energies, the fit does seem to match the spectrum 
relatively well. There is a tentative suggestion that the 
actual data is steeper than the fitted power law, but it 
appears to be a minor effect. 

We attempted to investigate this issue further by per- 
forming power law fits to the hard X-ray data only, 
thereby circumventing any problems introduced by the 
multitemperature disk blackbody model. We found that 
the power law index did appear to systematically increase 
by a value of ~ 0.3 to 1 when this procedure was used, 
but we were unable to obtain high enough signal-to-noise 
to reliably investigate the variability properties. There- 
fore, while we can conclusively state that the high energy 
power law component changes much more dramatically 
in 1997 than 2002, our data showing that this change 



specifically occurs in the power law index are less robust. 

Finally, although the power law is much more variable 
in 1997 than 2002, there does appear to be a small peak 
in the fitted power law index during the first oscillation in 
2002 (Figure[7|D) . There is also tentative evidence for this 
behavior in the first oscillation of other X-ray cycles in 
our 2002 observations. These detections are very close to 
our signal-to- noise limit, and thus we cannot determine 
whether the power law component varies in all the 2002 
oscillations but is in general too weak for us to detect 
(since the subsequent oscillations are not as strong as the 
first one), or whether the first X-ray oscillation is in fact 
physically distinct, and produces power law variations 
even though the others do not. 

6. ORIGIN OF THE "CLASS C" SUBFLARES 

It is interesting to speculate on the origin of the ~ 5-10 
mjy "class C" infrared flares which seem so ubiquitous 
i n GRS 1915-H05. 

lEikenberrv et al.1 1)200(1 ) showed that when these fiares 
were observed in July 1998, each flare preceded its asso- 
ciated X-ray oscillation by ^ 200-600 seconds, including 
examples in which the infrared light curve peaked and 
started to decay before the X-ray disruption began. To 
see if this property applies to the events we observed in 
2002, we show in Figure El a closer view of the isolated 
"subflare" from our first night of observations. The X- 
ray oscillation associated with this subfiare consists of 
a smooth rise in the accretion disk fiux, followed by a 
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sharp spike. As can be seen from Figure O the subflare 
peaks after the X-ray spike, but its initial rise does seem 
to begin ~ 30 seconds before the spike. This evidence 
is tenuous, but it is another hint that the ^ 5-10 mJy 
flares in GRS 1915-1-105 might all have a similar cause. 

Note that the X-ray event shown in Figure El is in no 
way unusual except for the long delay between it and the 
previous oscillations; a similar spike is seen at the end of 
most, if not all, of the 2002 oscillations (for example, see 
the last event in Figure Et>) ■ Similar oscillations have 
been discussed bv Taam. Chen fc Swank l,1997j|) . where 
the spike is referred to as a "secondary burst." Spectral 
fitting of these spikes is inconclusive but indicates that 
they may consist of a momentary increase in the flux from 
the accretion disk and power law components of the spec- 
trum, which makes them different from any other part 
of the 1997 or 2002 oscillations. Not e that the "class 
C" inf rared flares originally defined bv lEikenberrv et al.l 
l)2000D were associated with "soft dips" in the X-ray light 
curve that occur in the 1997 oscillations (see tl5.3|l . but 
for the purposes of this paper, we will refer to any in- 
frared subflare that occurs during the ~ 30 minute events 
and is associated with an individual oscillation as a "class 
C" flare. 

We cannot directly rule out a thermal origin for the 
2002 subflare (since it is an isolated event with no radio 
coverage), but the near-simultaneous presence of radio 
flares during our observations (see i|3.2|l . combined with 
the idea that all of the infrared flaring behavior in 2002 
was due to superimposed smaller events (see fQ, strongly 
suggests that each small infrared flare has a radio coun- 
terpart, and theref ore a nonthermal origin ( see, however, 
Jflljl . In addition. lEikenberrv et all l|200nli argued that 
the July 1998 "class C" infrared flares were probably re- 
lated to synchrotron jet emission, based on observations 
by Fcroci ct al. ( 1999) of a delayed, ^ 40 mJy radio coun- 
terpart to an X-ray oscillation spectrally similar to the 
July 1998 oscillations (again, however, we point out that 
the 2002 oscillations do not share this similarity). 

It is clear that the "class C" flares are a complex phe- 
nomenon, with a variety of X-ray counterparts. There is 
also further fast variability within the "class B" events 
that has not been discussed here: large (up to ~ 70 mJy) 
subflares in some of the 1997 data that unfortunately 
have no X-ray coverage, and smaller 10-20 mJy) 
flares that do not occur during the X-ray oscillations but 
rather during the initial rise of the infra red light curve, 
right after the trigger spike (Rothstein fc Eikenberrvl 
I2OOI . If we accept a synchrotron jet origin for all the 
flares in GRS 1915-1-105, then the observations discussed 
here show that the accretion disk can interact with the 
jet under a wide variety of timescales and circumstances. 

What is the nature of this interaction? The fact that 
at least some of the flares precede their associated X- 
ray oscillation argues for an "outside-in" origin for these 
events, in which the flare occurs somewhere outside of 
the inner accretion disk and eventually propagates in- 
ward to cause a disturb ance in the X-rays emitted there 
ijEikenberrv et alJl2000D. Simulations of magneti z ed ac- 
cretion disks bv lLovelace. Romanova fc NewmanI l|1994 fl 
have shown that an increase in the poloidal magnetic fleld 
in the outer part of the disk can lead to an enhanced ac- 
cretion rate, which triggers a slightly delayed burst of 
radiation from the inner disk. An alternate idea is that 



the flare could be produced by a reverse shock originating 
within a continuous jet outfl ow, which then propag ates 
back down to the inner disk ijEikenberrv et al1l2000|) . 

The possibility that shocks within a preexisting 
jet could be responsible for the infrared flares from 
GRS 1915-1-105 is intriguing in light of the variabil- 
ity that we find within each flare. A model for in- 
ternal sh ocks in m icroquasar jets was presented by 
Raiser, S unvaev fc S pruit (2000). It was primarily used 
to explain the larger "class A" events, but a modi- 
fied version might be applic able to the ~ 30 minute 
"class B" events as well (see iTiirler et all I2004j) , espe- 
cially given that a continuous jet has been observed 
in GRS 1915-1-105 and i s known to be associated with 
the ~ 30 minute cycles fPhawan. Mirabel fc Rodrigued 
|2000() . In addition.'Co Uins. Kaiser fc Co x (2003) showed 
that parti cularly large (^ 300 mJy ) "class B" events ob- 
served by iFender fc PoolevI 1)20001) cannot be explained 
by a simple model in which the flare results from an in- 
crease in the injected particle density at the base of the 
jet. 

7. DISCUSSION 

In this section, we outline our scenario for the ^ 30 
minute variability cycles in GRS 1915-1-105, based on the 
observations presented in this paper. The key feature of 
this scenario is that although GRS 1915-1-105 always dis- 
sipates energy in the jet (via infrared and radio flares) 
when the accretion disk becomes active, this emission is 
at least partially due to a superposition of smaller events, 
one for each oscillation in the X-ray light curve. Some 
of the events in GRS 1915-1-105 that would previously 
have been described as "class B" flares (in the sense that 
there is one large ejection every 30 minutes) may re- 
ally just be a superposition of smaller "class C" flares, 
which blend together and mimic a larger event. Episodes 
in which a large ejection really does occur, leading the 
infrared flare to dissociate from the X-ray light curve 
as the ejected material becomes causall y separated from 
the inner disk l|Eikenberrv et al.lll998a|) . seem to require 
a "trigger" spike in the X-ray light curve at the moment 
of the state transition, when GRS 1915-1-105 does not 
immediately begin the fast X-ray oscillations but rather 
enters a phase of duration ~ 300 seconds in which the 
accretion disk is hot and visible but changes on a much 
slower timescale. 

7.1. A Continuum of Behavior 

It is important to realize that the ~ 30 minute cycles in 
GRS 1915-1-105 do not neatly fall into discrete categories. 
From perusing the figures in this paper, one can easily 
imagine a continuum of behavior, from an extreme exam- 
ple in which the initial ejection completely dominates the 
infrared flare (Figure |3Ji) , to examples where the initial 
ejection is weaker and the excess emission from blended 
together faint flares is easier to see (Figure |3Jd) , to ex- 
amples from our 2002 observations where most of the 
infrared emission is due to blended together faint flares 
but a small "peak" is still observed at the beginning of 
the event, perhaps indicating a slightly larger initial ejec- 
tion (some of the flares in Figure|2li), to examples on the 
other extreme, where there is no peak at the beginning 
of the flare and the event is consistent with the initial 
ejection being exactly the same as the subsequent ones 
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(Figure int)- 

In this vein, we point out that X-ray observations 
which contain a strong trigger spike (the "class /3" state 
defined by Belloni ct al. 2000) have been observed to 
be accompanied by infrared flares with a wi de range of 
amplitu des. The events observed by Eikenb^rrv et all 
l)1998a(l clustered around ~ 100 mJy but ranged from 
^ 60 to 2 00 mJy, while a single event observed by 
IMirabel et al. (1998) during this X -ray state only reached 
a peak amplitude of 30 mJy. Meanwhile, the superim- 
posed faint flares reported in this paper during the "class 
q" state in July 2002 peak relatively consistently at ~ 30 
mJy, but even these less extreme events may have some 
va riation in ampl i tude; a single "class a" event observed 
bv lMirabel et alj l)1998|) was accompanied by a flare that 
only appeared to peak at 10 mJy. 

The X-ray behavior during the ~ 30 minute cycles may 
also represent a continuum, between cases where there is 
a clear X-ray trigger spike and cases where this distinc- 
tion is not so clear. Although the initial X-ray oscillation 
in each 2002 series is very similar to the subsequent os- 
cillations, it does share a few subtle characteristics with 
the 1997 trigger spike. For example. Figure ^ shows 
that the initial 2002 oscillation is "double-peaked," and 
at the end of the first peak the accretion disk flux drops 
sharply while the inner disk temperature continues its 
gradual rise. This is qualitatively the same behavior that 
was observed during the 1997 trigger spike (see and 
thus the first peak in this oscillation may be regarded 
as a weak analog of the trigger spike in 1997, albeit one 
which leads immediately to the X-ray oscillations rather 
than requiring a ~ 300 second phase in which the accre- 
tion disk is much less variable. Based on these proper- 
ties, it is not surprising that the initial 2002 oscillation 
sometimes appears to be accompanied by a slightly larger 
infrared fiare (Figurel^L, as discussed above) and slightly 
more variation in the power law index (see ti5.3|l than the 
oscillations which follow it; it may really be a weak vari- 
ant of the 1997 spike which leads to the large ejections. 

7.2. When Does the Ejection Occur? 

The observations prese nted in this pape r cont radict 
some of the conclusions of lKlein- Wolt et alJ l|200l . who 
analyzed simultaneous X-ray and radio observations of 
GRS 1915-1-105 and argued that the radio flares are pro- 
duced via a continuous ejection of material during the 
long, spectrally hard dips i n the X-ray light cu rve ("state 
C" in the classification of iBelloni et al.l[2000() . 

Independent of any of the specific conclusions we have 
reached in this paper, our observations pose a serious 
challenge for this argument. As is easily visible from Fig- 
ure[21 the duration of each infrared fiare in our July 2002 
observations is correlated with the duration of its accom- 
panying period of X-ray oscillations, not with the previ- 
ous X-ray dip. This strongly suggests that the fiares are 
associated with the X-ray oscillations, while providing no 
indication that they are infiuenced by any properties of 
the preceding d ip, as might be expected in the model of 
IKlein-Wolt et a l. ( 2002) , where a longer dip means that 
more material is being ejected into the jet. 

The differences which we observe between the 1997 
a nd 2002 infrared fiares also argue against the model 
of lKlein- Wolt et alJ l|200^ . The X-ray dips in 1997 and 
2002 are very similar to each other, and they are sepa- 



rated by X-ray oscillation cycles which have similar du- 
rations. The only real difference is that the 2002 dips 
appear to be slightly longer (~ 1000 seconds as opposed 
to ~ 600 seconds) and to have disk blackbody emission 
that is weaker by a factor of ^ 2 (although it is impor- 
tant to note that the disk blackbody emission is so weak 
in both cases that it is barely detectable) . It is not clear 
how these slight differences in the X-ray dips could lead 
to dramatic diffe r ences in the infrared flares under the 
IKlein-Wolt et alJ l|2002() model. Therefore, the fact that 
there appears to have been a large ejection during each 

30 minute cycle in 1997 but not in 2002 argues for a 
scenario in which the flares are associated with proper- 
ties of the X-ray state transition (such as the presence 
or absence of the trigger spike) and subsequent accretion 
disk oscillations, which do proceed much differently in 
1997 than in 2002. 

We do not dispute the conclusion of IKlein-Wolt et all 
(pOO^ that the presence of "state C" dips is required 
for infrared and radio flares to occur; for example, it 
is definitely possible that mass which is later avail- 
able to be ejected steadily builds u p during the dips 
(e.g. IBelloni. Midiari fc Fended l2000|) . It is also possi- 
ble that matter is being continuously ejected during the 
dips but that this material does not directly produced 
the large flares; instead, the flares could be produced 
via internal shocks which occur in the continuous jet 
when faster-moving material is ejected at the moment 
of the state transi tion (Vadaw ale et al.l|2003t see also 
[Fender. Belloni fc Gallo 2004^) . The main conclusion of 
IKlein-Wolt et al. (2002) which we refute is simply that 
the fiares can be entirely explained by the evolution of 
particles which were ejected during the dip, without ap- 
pealing to some other event that occurs afterwards in the 
accretion disk. 

The above ideas were hinted at by previous 
multiw avelength obser v ations containing i n frared 
data llEikenberrv et af] Il998at IMirabel et all 119981: 
lEikenberrv et al.ll2000|) . but they are confirmed by our 
present work. Nonetheless, it will be important to test 
them with future radio observations. The observations 
analyzed bv lKlein-Wolt et alJ (|2002j) did not contain any 
examples of the "class a" X-ray state for which there 
was significant overlap in the X-ray and radio coverage, 
and in general, there are only a few radio observations 
in the literature which overlap with the "class a" state. 

7.3. Infrared and Radio Flares 

One of the few previous multiwavelength observations 
during the "class a" state raises some questions about 
our current results. Ueda et al. ( 2002) obtained ~ 3 
hours of X-ray and radio (and, separately. X-ray and 
infrared) observations during this state. They detected 
at least one radio fiare that had a shorter duration than 
its accompanying X-ray oscillation period, with the main 
part of the flare decaying several hundred seconds before 
the end of the oscillations. This flare was most clearly 
seen at 1.3 cm (23 GHz), where i t reached an amplitude 
of - 15 mJy (Figure 2 of lUeda et al. 2002). 

In most models, radio flares are expected to 
peak later and last longer than their associated in- 
frared flares, due to opt ical depth effec ts in the 
ejected material (Ivan der Laan 1966: Mirabel et alJl998t 
IGollins. Kaiser fc Goxll20'o^ Therefore, the radio flare 
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observed bv lUeda et al.1 1)2002() does not appear to be the 
counterpart of the superimposed infrared flares we r eport 
in our current work, and, in fact. lUeda et al.l l)2002D used 
their observations to specifically rule out the possibility 
that the radio flares are produced via a superposition of 
smaller flares associated with each X-ray oscillation. In 
addition, the baseline flux density of the radio emission 
during these observations was 10 mJy; thus, if there 
was any "radio excess" associated with the X-ray oscil- 
lations after the decay of the flare, it was limited to, at 
most, this value. 

How can we understand these results? If the 
iTJeda et al.l l)2002 l) observation is representative of the ra- 
dio behavior during the "class a" state, then the initial 
ejections at the moment of the state transition may really 
be distinct, even though they do not appear that way in 
the infrared light curve. The initial ejection might have 
a relatively flat infrared-to-radio spectrum (and thus a 
prominent radio signature), while the repeating events 
associated with each X-ray oscillation could have little 
associated radio emission. If this is the case, then the 
differences between the "class a" and "class /3" states 
are not really as strong as we have indicated in this pa- 
per; both might produce ejections at the moment of the 
state transition which are substantially different from the 
subsequent ones (although it still the case that the "class 
/3" events produce, on average, much larger infrared and 
radio flares than the "class a" events ) . 

Alternatively, it is possible that the lUeda et alJ l|20fl2f) 
event was not representative of the normal behavior dur- 
ing the "class a" state. We have analyzed the X-ray light 
curve and spectrum from this event and found that there 
is a relatively sharp double peak in the flrst oscillation 
in this series, whose flrst peak has spectral and tempo- 
ral characteristics particularly reminiscent of a "class /3" 
trigger spike (see H7.1|l . In addition, there is an unusual, 
~ 10 second spike in the X-ray light curve which occurs 
~ 30 seconds before the state transition (although it is 
spectrally different from the trigger spike and appears 
to simply be a thermal "bimrp" in the spectrum). This 
behavior suggests that Ucda ct al. (2002) may have ob- 
served an anomalous event, in which the initial ejection 
was much different than those which normally occur dur- 
ing the "class a" state. 

Furthermore, the only other "class a" radio observa- 
tion we could flnd in the literature appears to be con- 
sistent with the in terpretation we have put forth in this 
paper. Figure 1 of iMirabel et all l|1998f) shows that dur- 
ing the "class a" state, the flare duration increases with 
wavelength as expected, including one case in which the 
flare was simultaneously observed in the infrared and 
three different radio wavelengths. The X-ray oscillation 
period associated with this event was not observed, but it 
is constrained to have been within a short (~ 30 minute) 
gap in the X-ray coverage that overlaps with the infrared 
flare and is of signiflcantly shorter duration than the ra- 
dio flares. Thus, it is broadly consistent with our current 
observations. 

The ultimate answer to the questions raised here may 
rest on how much variation there is in the amplitude 
of the "class C" flares associated with eac h X-ray os- 
cillatio n. If the "radio excess" seen in the lUeda et al.l 
l)2002tl observations (< 10 mJy if it exists at all) is typ- 
ical, then it is difflcult to see how the repeating, larger 



ra dio flares during the "class a" state (see H3.2I as well 
as IMirabel et alJfl998D . could be explained without re- 
sorting to initial ejections that are different from the 
subsequent events. Also note that the "class p" state 
in GRS 1915+105, which consists of an extremely long 
series of oscillations very similar to tho se seen during 
the "class a" state l|Belloni et al.l l2000() , is associated 
with stea dy radio emission of only ~ 4 mJy at 15 GHz 
l|Klein-Wol teral] [20021 . 

To further complicate the picture, the infrared obser- 
vations of Ucda ct al. (2002) show at least one instance 
(and likely more) in which a delayed ^ 30 minute infrared 
flare is seen during the "class a" state, beginning at the 
end of its associated X-ray oscillation period rather than 
the beginning. (Note that as in our observations, the 
duration of this flare still matches the duration of the 
oscillations.) It is unclear what could have caused this 
unusual delay, which is inconsist ent with a ll the infrared 
flares we report in this paper. lUeda et "a l. (2002) sug- 
gested the possibility of an internal shock occurring far 
out along the jet, but it is a puzzling event under any 
interpretation. 

Clearly, the only way to resolve the issues raised here is 
to obtain simultaneous X-ray, infrared and radio observa- 
tions over a long stretch of events. In the absence of such 
observations, we can only point out that the conclusions 
we have reached in this paper apply to all ~ 15 simulta- 
neous episodes that we observed. Overall, it is clear that 
the 1997 "class /9" events produce qualitatively differ- 
ent flaring behavior than the 2002 "class a" events and 
that these differences are related to the different ways 
in which the state transition proceeds and in particular 
to the X-ray trigger spike, which appears most strongly 
in the "class /3" observations. However, the exact man- 
ner in which properties of the trigger spike control the 
behavior of the initial flare — for example, whether the 
"sharpness" of the spike or the delay between it and the 
subsequent oscillations is the most important feature — 
cannot yet be conclusively determined. 

8. CONCLUSIONS 

We have presented simultaneous infrared and X-ray 
observations of the microquasar GRS 1915+105 (at 1- 
second time resolution) during a period of jet ejection in 
July 2002. By comparing these observations to those ob- 
tained in August 1997 by Eikenbcrrv ct al. (1998a|), we 
arrive at a picture in which a large ejection may not al- 
ways be present during the ^ 30 minute state transitions 
in GRS 1915+105. The observed infrared flares may 
sometimes entirely be a superposition of smaller, more 
complex phenomena, with large, infrared-bright ejections 
superimposed on them only when there is a "trigger" 
spike in the X-ray light curve. In particular, we find 
that: 

1. The duration of each infrared flare in 2002 matches 
the duration of its accompanying X-ray oscillation 
period. 

2. There is one instance in which a single, isolated 
X-ray oscillation occurs in the light curve, accom- 
panied by a faint infrared "subflare" superimposed 
on one of the main flares. 

3. These observations are consistent with a scenario 
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in which each X-ray osciUation has an associated 
faint infrared flare and these flares blend together 
to form, and entirely comprise, the larger events. 

4. The main difference between our X-ray observa- 
tions and the August 1997 observations is the pres- 
ence of a sharp "trigger" spike in the 1997 X-ray 
light curve during each ~ 30 minute cycle, which 
occurs ~ 300 seconds before the accretion disk os- 
cillations and signals a sharp change in the spectral 
properties of the source. 

5. The trigger spike appears to be associated with 
initial, large infrared flares, which are observed in 
1997 but not in 2002. 

6. The 1997 X-ray oscillations contain dramatic 
changes in the high energy power law component of 
the X-ray spectrum, which are not observed during 
the 2002 oscillations. This is somewhat puzzling, 
since both sets of oscillations produce similar in- 
frared subflares. 



IYadavll200lt iFender. Belloni & rTallcll200l . the obser- 
vations presented in this paper conflrm it. Our 2002 
X-ray obs ervations, whic h have many similarities to the 
[Eikenberr y et all l)1998afl observations except for the ab- 
sence of a distinctive trigger spike, do not show the same 
large flares in the infrared. Thus, these data are impor- 
tant in allowing us to pinpoint which aspects of the X-ray 
behavior are related to which types of events in the jet. 

Overall, we flnd that jet production in GRS 1915-1-105 
is a complex, unsteady phenomenon, down to the fastest 
timescales that have been observed. The jet and the 
accretion disk appear to interact on a regular basis, per- 
haps every time an oscillation occurs within the inner 
portion of the disk. These details need to be taken 
into account in accretion disk models, such as in the 
simulati ons of radiation pressure dominated dis k evolu- 
tion by iJaniuk. Czernv fc Siemiginowskal l)2002() , which 
found that energy dissipation in a corona or via an out- 
flow from the disk may be important in reproducing the 
GRS 1915-fl05 X-ray emission. 



7. There is tentative evidence that the infrared sub- 
flare observed in 2002 leads its associated X-ray 
oscillation by ^ 30 seconds, possibly suggesting an 
"outside -in" origin for these events, as in observa- 
tions bv lEikenberrv et alJ (j2003)- 

8. It is unlikely that any of the 2002 flares can be 
explained solely by the continuous ejection of ma- 
terial during the spectrally hard X-ray "dip" which 
precedes each set of X-ray oscillations, as suggested 
by Klein- Wolt et al. (2002). 

Although many authors have speculated on the impor- 
tance of th e trigger spike as a signature of large iet eiec- 
tions (e.g. lEikenberrv et al.lll998at iMirabel et all 119981: 
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